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In previous publications ||T]-§J, which have generated more attention than an- 
ticipated f§-fl, we showed that the n-fragment multifragmentation probabili- 
ties at any given transverse energy P n (E t ), could be reduced by means of the 
binomial expression 

O 



^ ■ " n\(m — n)\ 

to a single one fragment probability p(E t ), where E t = J] E; L sin 2 9i, the sum 
of the kinetic energies of all charged particles in an event weighted by the sine 
squared of their polar angle. This "reducibility" with E t was taken as a strong 
^ ■ indication that fragments were emitted nearly independently of one another 

with the same probability p(E t ). We further noticed that, under the assump- 
tion that Et is proportional to the excitation energy E (and thus temperature 
T oc y/Et), the Arrhenius plots of log(l/p) versus 1/y/Ei, produced straight 
lines. We called this "thermal scaling" because it suggests the Boltzmann form 
for the elementary probability p oc exp(— B/T) or 

, IBB' 
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P T y/Wt 



1 The Uroboros dragon, or the snake eating its own tail, is an archetype that can 
be traced from the Egyptians through the Greeks and the Gnostic religions into 
the Alchemical writings of the Middle Ages and Renaissance all the way to modern 
times. A symbol of self reference (think of Escher's hand drawing the hand that 
draws it...), it has been used to explain the universe. It has been applied by the 
authors of the preprint with a somewhat more modest intent. 
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with B representing the barrier for fragment emission. 



The present preprint attacks both of our claims. The first on the basis that our 
binomial fits could be matched in quality by a constrained Poissonian function. 
The second on the basis that the Arrhenius plots are allegedly shown to be 
the result of a strong autocorrelation. 

Since the second criticism is potentially more devastating than the first, we 
begin with it. 

As a preemptive strike, the authors raise concerns about the relevance of using 
\[Ei as a measure of the temperature. They cite as evidence to the contrary 
recent work exploring the impact parameter dependence of temperatures ex- 
tracted from excited state populations fn|. In that work a value of T=4.6 
MeV is determined, largely independent of impact parameter for moderately 
central to very central collisions of Au+Au at 35 MeV/nucleon. We are not 
unduly concerned by this result since temperatures of 4-5 MeV are almost al- 



ways extracted from such studies, independent of bombarding energy JTT|,[T2[, 
entrance channel mass asymmetry [0,0 > as wen as impact parameter fT0| , [I~3"| . 
Since such extracted temperatures seem to be insensitive to most features of 
heavy-ion reactions, we would not be surprised to find that they are also in- 
sensitive to E t . 

Furthermore, temperatures extracted from the double isotope ratio method 
I3| do indeed show a broad range of temperatures for heavy ion reactions at 



intermediate energies |T5 . 



Coming to the heart of their criticism, the authors argue that the observed 
linearity of the Arrhenius plots is due to autocorrelation. However, a fa- 
tal flaw mars the derivation of the proposed explanation for linear Arrhe- 
nius plots H] and undermines the conclusions of the preprint. This flaw is 
circular logic in the algebra, leading to a trivial identity that always holds, 
regardless of whether or not there exists an autocorrelation. 

Assuming with the authors that the binomial number of "throws" m is con- 
stant (for many experiments this is a good assumption) we consider the Ar- 
rhenius plot constructed with the average intermediate mass fragment (IMF) 
multiplicity (n) = mp and E t . So we plot logl/ (n) versus 1/y/El and in all 
experimental cases we obtain a straight line. As we have observed ||, the au- 
thors also observe that E t is constructed both with the light charged particles 
and the IMFs. In the authors' notation 

E t = {N c -(n))E™ + {n)E?™, (3) 



where Nc is the number of charged particles in an event and E\" and El are 
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Fig. 1. Top panel: The transverse energy E% as a function of the average IMF 
multiplicity (n) for the reaction Ar+Au, E/A=110 MeV. The line is a straight-line 
fit, Eq. (H|). Bottom panel: (n) as a function of E t . The fit is the exponential fit of 
Eq. (|) used in refs. (^g. 

the average transverse energies of the light particles and IMFs, respectively. 
The use of (n) to construct this independent variable E t creates the danger 
of autocorrelation in the study of the dependence of (n) on E t (as in the 
Arrhenius plot). This we also have discussed extensively ||. 

Now the authors |9[] observe that (n) and Nc are globally very strongly cor- 
related. They fit this correlation with the linear form 

N c = a' + V (n) (4) 

and finally they express E t as a function of the single variable (n) 

Et = b {l + {n) ) (5) 
where a = a'E^ and b = (V - l)£ t LP + Ej MF . 

A quicker but nearly equivalent procedure would have been to fit the equally 
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strongly correlated quantities E t , (n) directly with the above form (see the 
top panel of Fig. |l|) since Nq and E t are very tightly correlated as well. The 
values of a and b from a direct fit using Eq. (|5]) are 100 MeV and 220 MeV 
respectively. These values are nearly equal to those extracted by the authors 
(a=106 MeV, 6=213 MeV) when starting with Eq. (g). 

A moment's reflection shows that Eq. (|5]) contains the original experimental 
dependence of (n) on E t (bottom panel of Fig. [I]), and it does not differ in any 
significant way, other than in the linear fit, from the data used to construct 
the Arrhenius plots. It is just its inverse! Here is the start of the vicious circle . 

In fact, to fit the data they could have used with better results the Boltzmann- 
like functional form which gives linear Arrhenius plots as we have done in the 
bottom panel of Fig. |1| One can argue on which of the two is a better fit, but 
not on the fact that, for the two functions, one is approximately the inverse 
of the other. 

Armed with Eq. (|5]), the authors take an experimental Arrhenius plot which 
is linear, and they substitute E t with the value calculated from Eq. (|5|). 
Since the linear form is a good approximation to E t , nothing much happens 
to the Arrhenius plot which of course remains linear . It simply means that 
the form of Eq. (|5|) is a good approximation to the data! 

Yet the authors make the surprising claim that the linearity in such a plot is 
the result of an autocorrelation because now one has (n) both in the ordinate 
(as logl/ (n)) and in the abscissa (as 1/ \jb{ajb + (n))). 

In what follows we show that this technique for demonstrating autocorrelations 
is so general and fundamentally flawed, that it will give the same result for 
all functions, even for those in which autocorrelation has been ruled out by 
construction. We call this method the Uroboros method, for reasons that will 
become clear below. 

As we noticed above, Eq. ([|), either fitted directly to the data or indirectly 
through Eq. ([|), gives an empirical approximation to the total dependence of 
E t upon (n) 

E t = 9({n)). (6) 

On the other hand, we have the total dependence of (n) upon E t directly from 
the data 

<n> = / {Et) . (7) 
Clearly g ((n)) is an approximation to the inverse of / {E t ), g f~ x . Now we 
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construct the Arrhenius plot 



F((n)) = F(f (E t )) (8) 

with 

(n) = f (E t ) oc exp ^) . (9) 

If the Arrhenius plot is linear it means the f'(E t ) = f. Hence, 

F «n» = F (/ (Et)) =F(f(g «n»)) = F (/ «n)))) = F «n» (10) 

which is an identity, in the limit of g = f~ l . 

The vicious circle of circular logic can be summarized as follows: given any 
function y = y(x), substitute in it x with its value x(y): y = y(x(y)). This 
gives always the identity y = y. 

This result proves, if there was any doubt, that any function is completely 
correlated with its own inverse, and it is a great source of straight lines, irre- 
spective of whether there is an autocorrelation or not. 

The flaw, in our specific case, lies with the fact that Eq. @ gives the total 
empirical dependence of E t on (n) irrespective of whether there is autocorre- 
lation or not and consequently it is the approximate inverse of Eq. (|7p and of 
Eq. (j^). Hence the Uroboros, the snake eating its tail. 

It is left for the diligent reader to test the universality of this method by 
applying the Uroboros to discover autocorrelations in the most venerable de- 
pendences in physics. 

It is clear from the previous considerations that the authors' vicious circle 
does not prove anything, least of all autocorrelation. 

In order to appreciate the true role of autocorrelations, let as before 

(n) = f(E t ) (11) 



and 

Et = Et (E, (n)) (12) 
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where E is the excitation energy of the system. The dependence of (n) on E t 
can be best explored differentially 



a / \ d in) . ^ 
A in) = -±J-AE t . 
W dE t 



(13) 



Now, the change in transverse energy can be written 



AE t 



dE t 



BE 



AE + 



dE t 



d (n) 



A(n) 



(14) 



where the first term represents the change in E t due to the change in exci- 
tation energy and the second term represents the trivial change in E t due 
to autocorrelation. The partial derivative JP| at constant E (not the total 
derivative) tells us about the autocorrelation. 



Substituting Eq. ( |T4"D in Eq. ( |T3D we have 

1 d(n) dE t 



A (n) 1 



dE t d(n) 



d(n) dE t 



dE, dE 



AE 



(15) 



(n) 



where the complete left hand side is the increase in (n) induced by AE 
corrected for the autocorrelation . Failure to appreciate the difference between 



dEt 
d(n) 



and leads to disastrous circular reasoning. 



In fact, if we write 



AE t 



dE t 
din) 



A(n) 



(16) 



and substitute in Eq. ( 13|) we obtain as the authors did 



. . . din) dE t . . . . , , 
A In) = }J , . . A in) = A in) 
w dE t d(n) x 1 w 



(17) 



or the Uroboros (again). 

Of course, the authors might argue that they know the first term of Eq. ([14]) 
to be close to zero because they believe the temperature changes at most 
by 1 MeV. Then indeed only the second term, or the autocorrelation would 
survive. This position is clearly untenable. In fact, were the authors correct, 
the derivative in Eq. (|16"D (which would now be the same as that in the second 
term of Eq. (|T4])) should have the value of approximately 20-40 MeV per 
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IMF, which is the measured mean transverse energy per IMF fTq , [T7|l for this 
particular system. Instead, it has the enormous value of 220 MeV per IMF! 
(See Fig. |I|.) This alone should make the authors wonder and worry about 
constant temperature. In any case, as we pointed out above, their method will 
always give the same result, independent of the presence or abscence of an 
autocorrelation. 

For those who like to touch with their hands, here is a concrete example 
from an event generator code || which, by construction, generates complex 
fragment events according to 

p = exp (-B/T) (18) 



with 



T 



(19) 



where T is the temperature and E is the excitation energy of the system. 
In other words, the linear Arrhenius plot is contained in the events by con- 
struction. By using the independent variable E instead of E t , we are assured 
of avoiding the autocorrelation problem altogether. 

We construct the Arrhenius plot 1 / (n) versus 1 / \fE which turns out to be 
linear as expected (circles in the bottom panel of Fig. Q). Then we notice that 
E and (n) can be plotted one against the other (circles in the top panel of 
Fig. 0) and we approximate 

E = a + b (n) (20) 



by fitting the dependence with a straight line. Then we substitute E with the 
value given by Eq. ( |211D and note, as the authors, that the Arrhenius plot con- 
structed from this substitution (solid line in bottom panel of Fig. |2|) remains 
linear over the range in E t where Eq. ( |2D| ) well represents the simulation. 
Hopefully by now the game has become rather tedious, and nobody would 
argue that the result is just due to an autocorrelation. 

To verify the extent of autocorrelations in this model, we follow the same 
procedure as above but use E t , like in the data. Now, we expect some auto- 
correlation. The question is: how much? We obtain a linear Arrhenius plot, 
whose slope should coincide with the slope of the previous plot after scaling E t 
by 3/2. (The model assumes E t = 2/3E.) Indeed it does so to within better 
than 10 percent |J (open squares of Fig. |) This difference in slope is the 
extent of the autocorrelation problem. By redefining (n) through an increase 



7 





io" 1 1 1 1 1 1 1 1 

.00 .01 .02 .03 .04 .05 

E-^MeV- 1 / 2 ) 

Fig. 2. Top panel: The excitation energy E as a function of the average IMF mul- 
tiplicity from a binomial simulation |J of the decay of a Xe nucleus (circles). The 
squares show the same correlation but using for E a scaled value of E t , E = 3/2E t . 
Bottom panel: The Arrhenius plot 1/ (n) versus l/y/E (circles) and 1/ (n) versus 
l/y/3/2Et (squares). The solid line represents the (erroneous) "self-correlation" ex- 
planation proposed by the authors (see text). 

of the element threshold (Z th ) in the definition of IMF (IMF: Z tb < Z < 20), 
we can substantially reduce even this small autocorrelation. Yet the authors 
would conclude that in both cases, using E or E t , the linear Arrhenius plots 
are completely due to autocorrelation. 

An even better procedure, which we can apply directly to experimental data, 
is the following. We can redefine the IMFs to be the fragments of a single 
element, say Z=10. Now the contribution of the Z=10 fragments to E t is so 
small that it can be neglected and consequently autocorrelation is negligible 
as well. Because of the restriction of the definition of an IMF (one single Z), 
the probability of emitting an IMF becomes very small, and in fact, we observe 
that the distribution of multiplicities becomes Poissonian. The Arrhenius plots 
of log (n) versus 1 / \[E~ t remain linear and the linearity can be tested for each 
Z as in Fig. [3] without the worry of autocorrelation. The remarkable set of 
straight lines shown in Fig. [| together with many other sets obtained for 
all the experiments we have analyzed, shows that thermal scaling is a real 
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36 Ar+ 197 Au, E/A=110 MeV 
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Fig. 3. The average yield per event of the different indicated elements (symbols) as 
a function of 1/y/Et. The lines are linear fits to the data. 

physical feature quite safe from the dreaded fangs of the Uroboros snake. Yet 
the procedure employed by the authors will always yield an autocorrelation, 
regardless of whether one is present or not. 

We now return to the first criticism, on which there is little to say. 

The authors claim that there is no a priori reason for choosing a binomial 
instead of a constrained Poissonian, because they both fit the experimental 
multiplicity distributions equally well. We have no great argument with this. 
We actually like the Poissonian case and we would have chosen it if the data 
had not indicated a very high probability and a constraint (in size or tempo- 
ral) to boot [[J. A binomial satisfied both requirements. By requiring a fixed 
number of tosses m, it incorporates both a scale and a constraint. There are 
other ways to achieve the same goal. For instance, we could chop the Pois- 
sonian above a prescribed n or do something else, like the authors did. In 
point of fact, when we restrict the definition of complex fragments to a single 
atomic number Z (see Fig. |3|), the probability falls so much that the "binomi- 
als" become Poissonian on their own, removing us from the embarrassment of 
choice. 

But binomials, Poissonians, and constrained Poissonians, all share the feature 
of being constructed upon one fragment probabilities. All of them construct 
the n-fragment probability P n from the one fragment probability p. All are 
reducible and this is our experimental point! A Poissonian is quintessentially 
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Fig. 4. The empirical scaling of mZ t ^ for different Z t h (different symbols) for the 
reaction 36 Ar+ 197 Au at E/A=110 MeV. 

reducible. Because of the loss of scale, its elementary probability is arbitrarily 
normalizable and coincides with the mean. Remembering that the Poissonian 
originates from the binomial, we get (n) = mp without however, being able 
to disentangle m from p. On the other hand, if we introduce a constraint, in 
whichever expression, we get back the approximate scale m m Zq/ (Z) where 
Zq is the source size. 

An example of a mass constraint in our binomial fits is given by the approxi- 
mate scaling that we have observed [Q in the experiments: 



^z th -^th ~ constant 



(21) 



shown in Fig. f|. 

In view of the above, we do not see any merit to the authors criticism. Their 
function is reducible by construction, since it assumes a priori independent 
emission (see their appendix). Furthermore, their constraint reintroduces the 
scale, by the back door, that the Poissonian threw out the window. Reducibil- 
ity is our key statement, binomial or Poissonian, as the case may be. 

The authors seem to imply a big difference between the two, but the issue is 
reducibility. Clearly, our experimental claim that the n-fragment multiplicity 
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distributions are reducible to a single one fragment probability is not ever 
called into question by the authors. 

In summary, the present preprint || attempts to give alternative or trivial 
explanations for the reducible and "thermal" nature of the experimental IMF 
excitation functions reported in We have demonstrated that the pro- 

posed "self-correlation" explanation for linear Arrhenius plots is based upon a 
flawed autocorrelation analysis involving circular reasoning. Equivalently, the 
authors have used a total derivative when they should have taken a partial 
derivative. This is a fatal mistake which produces the "self-correlation" result, 
independent of whether or not there are autocorrelations truly present. 
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